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Controlling electron spins strongly coupled to magnetic and nuclear spins in solid state sys- 
tems is an important challenege in the field of spintronics and quantum computation. We show 
here that electron droplets with no net spin in semiconductor quantum dots strongly coupled with 
magnetic ion/nuclear spin systems break down at low temperature and form a non-trivial anti- 
ferromagnetic spatially ordered spin-texture of magneto-polarons. The spatially ordered combined 
electron-magnetic ion spin-texture, associated with spontaneous symmetry-breaking in the parity 
of electronic charge and spin densities and magnetization of magnetic ions, emerge from ab-initio 
density functional approach to the electronic system coupled with mean-field approximation for the 
magnetic/nuclear spin system. The predicted phase diagram determines the critical temperature as 
a function of coupling strength and identifies possible phases of the strongly coupled spin system. 
This prediction may arrest fluctuations in spin system and open the way to control, manipulate and 
prepare magnetic and nuclear spin ensembles in semiconductor nanostructures. 

PACS numbers: 81.07.Ta,75.75.Fk 



There is currently significant interest in developing 
quantum information storage and processing capabilities 
[1] using electron (e) and/or hole (h) spins strongly cou- 
pled to spins of either magnetic ions (MI) and/or nuclear 
spins (NS) in a number of solid state systems. This in- 
cludes GaAs based gated two-dimensional [2 and zero- 
dimensional systems [3]-[5], InAs self-assembled quan- 
tum dots [SHE], CdTe quantum dots [9H2T], nanocrys- 
tals [22[ |23] , NV centers in diamond [24] , phosphor im- 
purities in silicon [25] and carbon nanotubes [26 . 

In these systems electron spins play either a role of 
qubits or are used to control Mis or NSs. For e-spin 
qubits, much effort is directed to determine the role of de- 
coherence by NSs [H |5] . For electron interacting with MI 
in diluted magnetic semiconductors [271 EH] , the central 
spin problem has been understood in terms of magneto- 
polarons (MP) - a cloud of magnetization surrounding a 
localized e-spin [28]. On the other hand, the interaction 
of many electrons in a spin singlet state with MIs/NSs, 
e.g., in a metal [29] or a quantum dot [Ml |15] in- 
duces Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion among MIs/NSs. The question of whether spin tex- 
tures could form in a strongly coupled two-dimensional 
electron-NS subsystem has been addressed recently by 
Loss et at. [30 . Both scenarios of MPs and RKKY inter- 
actions can be realized in semiconductor quantum dots 
containing MIs/NSs by changing electron/hole concen- 
tration with a gate, through modulation doping [TH [15] 
or through deformation of quantum dot confining poten- 
tial [17]. Spin singlet droplets are already employed in 
initialization of coded qubits in lateral quantum dots [3| 
and as component of trions in optical manipulation of 
spin in self-assembled quantum dots [7l|8]. 

Currently, magnetic ordering in closed-shell QDs 



doped with Mn ions is under question. In particular, 
recent work by Oszwaldowski et al. [21] suggests that 
closed-shell QDs doped with Mn do not allow magnetic 
ordering. The authors came to such conclusion because 
the ground-state of a nonmagnetic two electron system 
is a spin singlet but neglected RKKY interaction previ- 
ously discussed in Ref. [14 . Hence a full understanding of 
strongly interacting electronic singlet state with MI/NS 
in QDs still appears to be missing. In this work we show 
that the closed shell magnetic QDs strongly coupled with 
magnetic ions do allow magnetic ordering. The pertur- 
bative RKKY interaction among Mn ions for the closed 
shell QDs, first discussed in Ref. ^4 , does not change 
the symmetry of electronic ground state. In the present 
work, our nonperturbative approach predicts novel bro- 
ken symmetry electronic ground state with spin- textures, 
a non-trivial form of magnetic ordering, followed by fer- 
romagnetic state of electrons and magnetic ions at very 
strong coupling. Because of the degeneracy of the spin- 
textures under continuous rotational transformation, the 
magnetization of Mn and the charge and spin densities 
of e/h are expected to drift over time. We use ab-initio 
density functional approach to the e-system and mean- 
field approximation for the MIs/NSs. We show that be- 
yond RKKY approximation [14 a spin singlet e-droplet 
corresponding to closed-shell quantum dots strongly cou- 
ples with MIs/NSs. The e-droplet breaks down and 
forms a non-trivial anti-ferromagnetic (AFM) spatially 
ordered spin- texture of MPs at low temperature. For 
very strongly coupled and weakly confined quantum dot 
system this transition is followed by a transition to a fer- 
romagnetic state. This prediction opens the way to con- 
trol, manipulate and prepare MI/NS ensembles in semi- 
conductor nano-structures. 
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We focus here on closed-shell QDs |3T] containing 
A/' = 2,6,12,... electrons in the presence of many spins of 
either NSs or Mis, e.g., mangan in CdTe [14 } [T5 l l?7l[28] . 
We approximate many spins by a continuous magneti- 
zation p£l [17] and study the strongly coupled e-MI sys- 
tem as a function of temperature T, electron number 
strength of e-MI exchange coupling Jsd^ MI density rim 
and strength of confining potential ujq. In particular, 
we find that for a given confining potential, number of 
electrons, MI/NS density and their coupling with elec- 
trons, there exists a critical temperature, T* for a two- 
dimensional nucleation and growth of inhomogeneous 
AFM spin-textures with broken symmetry. Below T* 
and above a critical e-MI exchange coupling strength, 
spin singlet droplet and a homogeneous magnetization 
density breaks down and molecular states of MPs associ- 
ated with individual e-spins form. The MPs correspond 
to the inhomogeneous magnetic field of Mis inducing an 
effective spatially varying potential localizing electrons 
with spin up in different positions from the electrons with 
spin down, significantly changing electronic spin distri- 
bution. This is to be compared with closed-shell QDs 
with large confinement potential and with only two Mis 
O |15]. Using exact diagonalization [14] of the inter- 
acting e-e and e-MI Hamiltonian it has been shown that 
the RKKY coupling, a second-order effective interaction 
between Mis mediated by two e's with opposite spins, de- 
scribes well the ground state with total magnetic moment 
Mz = (AFM ordering) and = 2M, ferromagnetic 
(FM) ordering, depending on the relative position of Mis 
in a QD, while maintaining the spin singlet electronic 
ground state with spin polarization = 0. Here M is 
the total spin of single MI/NS, e.g., M = 5/2 for Mn. A 
related problem has been recently studied where an ana- 
lytical variational form of the two electron wave-function 
coupled to a large number of Mis, neither two-body sin- 
glet nor triplet, called pseudo-singlet, was introduced to 
describe partial quantum correlations of the coupled spin 
singlet-MI system [2T] . 

We use ab-initio density functional to describe the 
droplet of e's in a parabolic quantum dot with closed 
electronic shells for electron numbers = 2, 6, 12, .. and 
mean-field approximation for the MIs/NSs system |T6] . 
We employ spin unrestricted local density approxima- 
tion (LSDA) for electrons where the many-body Hamil- 
tonian is replaced by the Kohn-Sham (KS) Hamiltonian, 
HkS' In LSDA, the self-consistent KS orbitals are calcu- 
lated for spin up and down independently, without any 
additional symmetrization of their spatial dependence. 
The electrons interact with a magnetic field produced by 
the Mis which in turn is determined self-consistently by 
the electron spin density. The electrons fill KS-orbitals 
according to Fermi statistics at finite temperature and 
the electronic correlations are taken into account via 
exchange-correlation (XC) energy functional. In this ap- 
proach, the self-consistent solutions form a large class 



of variational many-body wave- functions, among them 
the configuration with the lowest free-energy is found. 
We also decompose the planar and perpendicular com- 
ponents of the confining potential of a single QD as de- 
scribed in Refs. [16J [17] and expand the electronic wave- 
functions in terms of its planar, ipicr{r)^ and subband 
wave function (,{z), and project Hks onto the XY-plane 
by integrating out (,{z), assuming that only the lowest 
energy subband is filled. Hence KS orbitals, tlji(j{f)^ are 
calculated by diagonalizing H^s'^iai'^) = ^iai^iai^i in 
real-space. Here e^cr is the KS eigen-energies, and Hks = 
+ Va with KS effective potential denoted as = 
Vqd ^ Vh + ^xc ~ ^^sdif)' ^ is the Planck constant, 
m* is the e effective mass, — e is the e-charge and a = ±1 
denotes spin up (t) and down {\). Vqd is the planar con- 
fining potential of QD, Vh and Vxc Hartree and spin 
dependent exchange-correlation potentials and hsd{r) = 
Jem J dz\^{z)\^BM{Mb{r,z)/kBT). Here Bm{x) is the 
Brillouin function [27 , /c^ is the Boltzmann constant, 
and b{ri) = Jsd[n^{'^i) — ^i(ri)]/2 is the effective field 
seen by Mis, and n^{Ti) = J2i Here 
/(e) = l/{exp[(e - iJi)/kBT] + 1} is the Fermi-Dirac 
distribution function, ji is the chemical potential and 
Jem = Jsd^mM is the mean-field e-Mn or e-NS exchange 
coupling. Jsd is the local exchange coupling, FM for e- 
MI and AFM for e-NS or h-MI. Sign of Jsd affects the 
relative orientations of e/h spin polarization with respect 
to MIs/NSs. Finally rim denotes MI/NS average density. 

Our approximation in neglecting spatial distribution 
of Mis is substantiated by experiments on colloidal 
nanocrystals and self-assembled QDs [22l [2^ showing 
that inhomogeneity in the distribution of Mn ions per 
QD, corresponding up to ~ 5% doping in typical mag- 
netic semiconductors, can be neglected [27]. For nuclear 
spins in, e.g., GaAs, inhomogeneity due to different iso- 
topes is unlikely to play an important role. Typical Mn 
density is much lower than the NS density, however, be- 
cause e-Mn exchange coupling is much higher than e- 
NS exchange coupling, Jgm for both can be of the same 
order of magnitude. Therefore we only present numer- 
ical results for a system of interacting e/h and Mis. 
For strongly coupled e-MI system we perform numeri- 
cal calculation for (Cd,Mn)Te where = 5.29 nm, and 
Ry* = 12.8 meV are the effective Bohr radius and Ryd- 
berg, the sd exchange coupling is Jsd = 0.015 eV nm^, 
effective mass m* = 0.106, and e = 10.6 ^2]. We con- 
sider QD's with electronic shell spacing ujq in the range 
of 1 — ?>Ry* ^ width of 1 nm, and variable MI density rim- 
For definiteness, we focus on the example of Mn isoelec- 
tronic impurity in CdTe, with the 2;-component of Mj of 
impurity spin satisfying = — M, — M + 1, . . . , M and 
M = 5/2. The direct Mn-Mn AFM coupling is negligible 
in the range of MI densities considered in this study. 

The spatial dependence of the z-component of magne- 
tization Mz{r) of Mn ions, solutions of the self-consistent 
LSDA equations for the coupled e-MI system, are shown 
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FIG. 1: The spatial profile of magnetization Mz for a closed 
shell parabolic QD with uq = ^Ry* and Um = O.lnm"^. (a), 
(b), and (c) correspond to N = 2,6,12 at T = 0.5K and 
Jsd — 15 meV nm~^. Coordinates [x^y) are expressed in 
effective Bohr radius. Inset in Fig.lc shows a decagon molec- 
ular structure oi N — 12 closed-shell QD with a molecular 
state of 12 MPs, each shown by an arrow. The direction of 
arrows show the direction of Mz in each MP. 



in Fig. [TJa-c) for closed shell parabolic QD correspond- 
ing to cjo = '^Ry* ^ rim = 0.1 nm~^ and temperature 
T = 0.5 K for N = 2, Fig. [ija), = 6, Fig. [ijb) and 
N = 12, Fig. [ijc). Note that the total magnetization 
per unit area A, (M^) = ^ / (Pr{Mz{r)) = 0, indicating 
net AFM ordering of Mis. These states clearly resemble 
spin textures with broken rotational symmetry. We see 
two spatially separated magnetization clouds for = 2 
quantum dot, six for = 6, and twelve for N = 12. 
Hence we consider these states as molecules of MPs. The 
inset shows schematically e-spin density for twelve elec- 
trons with a decagon ring structure and two electrons in 
the middle in a form of a spin-corral, a circular symmetric 
spin texture (see Fig. [s]). The inhomogeneous spin den- 
sity gives way to a uniform solution in which = 5^=0 
for T > T* ^ 2 K independent of TV. Note that this 
finding is in contrast with open shell FM states that are 
stable up to T ^ 20K [16, 17 and the FM state proposed 
recently in Ref. [21] for closed-shell QDs. 

To investigate the origin and stability of AFM states 
presented in this work, we focus on a closed-shell QD 



with N = 2. Let us first consider the RKKY inter- 
action between two Mis in a QD [14^ with = 2 at 
positions Ri (X, 0) and R2 = {-X,0): J{X) = 
—7(2 — 5X^)e~^ is the e- mediated effective interaction 
between two Mis, where 7 = (J^^/tt/q)^ x (1/(16cjo)) 
and distance is measured in Iq = ^/h/{2m*uJo). Here 
J^^ = Jsd^/d and d is the thickness of QD in perpen- 
dicular direction [TT. We see that for X > y^2/5 the 
interaction is AFM. It is now possible to imagine that 
the RKKY interaction triggers a broken symmetry state 
where the exchange interaction localizes an electron with 
spin up on the left MI and spin down electron on the right 
MI, a two-atomic molecule of MP, seen in Fig. [ija). In 
Figs. [2] we show self-consistent solutions of LSDA equa- 
tions including the effective potential Va- and spatial pro- 
file of spin-dependent KS wave-functions for parabolic 
confining potential with level spacing uq = 2Ry* . In 
Fig. |2|a), T = 2K and the paramagnetic state with 
Mz = Pz = is the ground state. The spin-dependent 
effective potential Va- is shown with open-circles in which 
= V^. In this case we find that the self-consistent 
KS eigen-energies and eigen-states with opposite spins 
are identical, eo^ = eo^ and tpo^ = tpo^. In Fig. we 
present our results for T = 0.5i^. Interestingly, with de- 
creasing temperature a continuous displacement in spin- 
dependent wave- functions develops that spontaneously 
breaks the circular symmetry of the QD as discussed 
above. In this case we find two degenerate spatially sep- 
arated wavefunctions ?/^ot V^oi with eo^ = eo|. This 
is analogues to the nucleation of para-H2 molecule from 
para-He if two nuclei of He undergo a continuous frag- 
mentation. A two-dimensional profile of e-spin polariza- 
tion Pz = (n^ — ni)/2 shown in (c) is also a result of 
spontaneous spatial spin-separation of KS orbitals. Two 
asymmetric bumps observed in effective potentials and 
Vi are the result of spontaneous formation of AFM spin 
texture in M^, shown explicitly in Figjl] 

The spontaneous symmetry breaking and the forma- 
tion of MPs depends on the parameters of the system. 
For example, in the example studied here, the critical 
temperature depends on e-MI exchange coupling strength 
Jem = Jsdrim^i which is proportioual to exchange cou- 
pling, MI density and MI spin. Figure [S] shows the phase 
diagram oi N = 2 closed shell QD with ujq = 2i??/* and 
Um = 0.1 nm~^. The curves show critical temperature as 
a function of Jem- In the limit of low Jem, e.g., either for e 
and/or low density of Mis, a direct transition from AFM 
spin texture, shown as inset, to normal state is seen. For 
high Jem, e.g., for either h and/or high density of Mis, 
a transition through an intermediate state of spin-corral, 
shown as inset (see Fig. |3| is observed. This is a transi- 
tion from rotational symmetry breaking state, stable at 
low temperatures, to another type of spin texture with 
rotational symmetry, stable within < T < . For 
Jsd = 75 meV nm^, corresponding to h-MI exchange 
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FIG. 2: The spatial dependence of wavefunctions ipo^^ipoi, 
effective potentials and Pz in a QD with A/" = 2, cjq = '^Ry* , 
rim O.lnm-^ In (a) T > T*, = = and V'ot V^o^- 
The dashed lines show the QD parabolic confining potential, 
Vqd. The spin-dependent effective potential Va = VQD + Vi/ + 
Vxc ~^^sd/'^ is shown by dashed lines with open circles, with 
V]- = V^. In (b) T < T* . A pronounced separation of spin 
up and down wavefunction and two asymmetric bumps in 
and are shown. In (b) the bottom of potentials and Vi 
are shifted to zero. The horizontal-axis used in (a) and (b) is 
plotted along the diagonal-direction of Pz shown as principal 
axis by a dashed-arrow in (c), and denoted by A in (b). Planar 
coordinates {x,y) are expressed in effective Bohr radius (ao). 




(meV) 



FIG. 3: Phase diagram, temperature T vs e-MI coupling 
strength Jem of N = 2 closed-shell QD with cjq = 2Ry* . 



N of electrons in a quantum dot, MI density, exchange 
coupling strength, temperature and confining potential 
is predicted. The predicted ordered spin states may ar- 
rest fluctuations in the spin system and open the way 
to control, manipulate and prepare MI/NS ensembles in 
semiconductor nanostructures for quantum information 
processing and storage. 
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coupling, and rim = 0.1 nm~^ we find = 17 and 
T2* = 35 K. We emphasize that the total spin of elec- 
trons is zero for the range of parameters considered in 
this study. However, by increasing Jgm, above a critical 
value that is determined by the energy difference in the 
lowest unoccupied and highest occupied KS energy levels, 
a transition to spin-triplet takes place [16] [TTl [19^. For 
the typical quantum dots investigated here, the critical 
value of J*^ ~ 30 meV is much higher than the range of 
exchange couplings Jem responsible for the spin textures 
predicted here. 

To summarize, low temperature solutions of ab-initio 
density functional approach to electrons in closed-shell 
quantum dots strongly coupled with the magnetic and/or 
nuclear spins described in the mean-field approximation 
show the existence of a broken-symmetry spatially in- 
homogeneous state of geometrically ordered molecules of 
magneto-polarons. The phase diagram of topologically 
stable non-trivial AFM states determined by the number 
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